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Abstract 

It is known that a large, charged body immersed in a solution of multivalent 
counterions may undergo charge inversion as the counterions adsorb to its sur- 
face. We use the theory of charge inversion to examine the case of a deformable, 
porous macroion which may adsorb multivalent ions into its bulk to form a three- 
dimensional strongly-correlated liquid. This adsorption may lead to non-monotonic 
changes in the size of the macroion as multivalent ions are added to the solution. 
The macroion first shrinks as its bare charge is screened and then reswells as the 
adsorbed ions invert the sign of the net charge. We derive a value for the outward 
pressure experienced by such a macroion as a function of the ion concentration 
in solution. We find that for small deviations in the concentration of multivalent 
ions away from the neutral point (where the net charge of the body is zero), the 
swollen size grows parabolically with the logarithm of the ratio of multivalent ion 
concentration to the concentration at the neutral point. 
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1 Introduction 



The chemical and biological function of a macromolecule, and of a polyelec- 
trolyte in particular, often depends critically on its ability to be stretched or 
compressed. The full genomic DNA of E. Coli, for example, must be com- 
pacted from its contour length of 1.4 mm (4.2 million base pairs) to fit inside 
a nucleolar region about 1/xm across p]. Achieving this kind of compaction 
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often requires biological organisms to expend a great deal of metabolic energy, 
but in some cases a macroion may be induced to spontaneously collapse into a 
tightly-packed condensate solely through the addition of multivalent counteri- 
ons to the surrounding solution. The trivalent ion spermidine, for example, has 
been known for several decades to induce the formation of DNA condensates 
[2] . A number of similar examples can be readily found for polyelectrolytes in 
a variety of contexts: the collapse of polyacid brushes [3lHf5] and the forma- 
tion of DNA precipitates j6], to name a few. In each of these cases, swelling 
of the polyelectrolyte can be induced or suppressed by adding/removing ions 
from the surrounding solution, thereby altering the polyelectrolyte's function. 
These last two examples, polymer brushes and DNA precipitates, provide par- 
ticularly interesting behavior, and they will be discussed more fully over the 
course of this introduction. 

One mechanism by which a macroion may be induced to swell or to shrink is 
through correlation-induced adsorption of multivalent ions. It is known that 
when a large, charged body is immersed in a solution of multivalent oppositely- 
charged ions, these ions can adsorb to its surface and form a strongly corre- 
lated liquid \7\. Such correlation- induced adsorption leads to screening of the 
bare charge which is much stronger than predicted by conventional Poisson- 
Boltzmann solutions. For a certain critical value of the counterion concentra- 
tion, the bare charge of the body is totally neutralized by adsorption, and for 
even larger concentrations additional adsorption inverts the sign of the net 
charge. This effect, dubbed charge inversion, has been described extensively 
over the past decade (a review can be found in Ref. [8]), and the theoretical 
framework that has been developed is used here to describe the swelling of a 
charged macroion in ionic solution. 

In particular, we extend the idea of charge inversion to the case of a large, 
porous, highly-charged molecule — like a coiled or globular polyelectrolyte — 
where the "bare charge" exists not just on the surface but within the bulk 
of the molecular "ball". In this case ions in solution can penetrate into the 
volume and adsorb onto some particular segment of the molecule [9]. If the 
bare charges within the ball are sufficiently dense and a sufficiently high con- 
centration of multivalent ions is present in the solution, adsorbed ions form 
a three-dimensional strongly-correlated liquid within the volume of the ball 
and cause it to undergo charge inversion. Thus far the theory of charge inver- 
sion has primarily been used to describe surface phenomena; it our intention 
to apply it to a situation like the one described here, where counterions are 
adsorbed within the volume of a charged body and cause it to swell or shrink. 

A number of recent experiments have observed swelling in systems of con- 
densed polyelectrolytes. The authors of Refs. [3] and [5] have studied swelling 
in polyacid brushes as metal cations are added to the surrounding solution. 
When salts with divalent or trivalent cations were added, like Ca 2+ or Al 3+ , 
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the brush height was seen to decrease by as much as a factor of ten. This 
behavior is qualitatively similar to what is known to happen with the ad- 
dition of monovalent salt. However, the studies of Refs. [313] and [5] report 
non-monotonic changes in the brush height as certain multivalent salts were 
added. The brush height was seen to collapse abruptly as multivalent salt was 
added, reach some minimum condensed size, and then swell again with the ad- 
dition of more salt. The non-monotonic change of height may be attributable 
to correlation-induced charge inversion. At some concentration of multivalent 
ions, the polyacid brush is neutralized and it attains its minimum size. If the 
ion concentration deviates from this neutral point, however, the brush is either 
over- or under-charged and it experiences an outward (expansive) pressure re- 
sulting from electrostatic forces. 

Similar behavior was observed in studies of the structure of DNA precipitates 
as the multivalent salts spermine (4+) and spermidine (3+) were added. The 
authors of Ref. [6J performed measurements of the interhelical spacing in liq- 
uid crystalline DNA precipitates induced by spermine and spermidine. In their 
study, short segments of nucleosomal DNA (~ 146 bp) were placed in ionic 
solution, and at some critical concentration of multivalent salt, the DNA seg- 
ments were observed to aggregate into tightly-packed condensates (bundles) 
of parallel, hexagonally-arranged DNA "rods". This arrangement is depicted 
schematically in the inset of Fig. [2j As spermine or spermidine salt was added 
to the solution, the interaxial spacing between segments in the precipitate 
bundles was seen to achieve a minimum (~ 28 nm, greater than the "hard 
core" diameter of double- helical DNA), and then increase again as more salt 
was added. Such behavior stands in contrast to the monotonic decrease in 
size observed to result from the addition of monovalent NaCl ions. We sug- 
gest that in this experiment the multivalent ions form a strongly-correlated 
liquid within the bulk of the condensed DNA bundles, eventually leading to 
charge inversion and non-monotonic swelling behavior as the ion concentration 
is increased past the neutral point. 

In this paper we consider correlation-induced adsorption in systems of porous 
macroions, and we suggest a simple description of non-monotonic swelling 
behavior for macroions in the presence of multivalent salts. We propose that 
at some concentration of multivalent counterions, a porous, charged macroion 
achieves a minimum size resulting from its bare charge being neutralized. As 
the ion concentration deviates from this neutral point, an outward pressure 
results from electrostatic forces. The way in which this pressure, and therefore 
the swelling of the polyelectrolyte, depends on the properties of the solution 
is not well understood. It is the purpose of this paper to explore overcharge- 
induced swelling via a simple model of porous charged bodies in solution, 
and to determine the functional dependence of swollen size on multivalent ion 
concentration near the neutral point. 
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Fig. 1. Schematic of the system to be studied. A porous, spherical, negatively 
charged body sits in a solution of positive multivalent ions (denoted by +Z's) and 
dissolved monovalent salt (+'s and -'s). Monovalent salt passes freely through the 
sphere while the Z-ions may adsorb to its volume, forming a strongly-correlated 
liquid. 

The following section details the physical model to be studied and presents a 
full solution for the pressure and swollen size based on a free energy description 
of the system. The Discussion section compares our theoretical predictions to 
existing experimental data and reviews the limitations of our model, followed 
by some concluding remarks. 



2 A simple model of swelling 



2.1 Model and approach 



We begin by considering a three-dimensional body in a solution of multivalent 
counterions and some monovalent salt, such as NaCl. We assume the body 
to be a sphere of radius R with some fixed charge Qq, which for the sake 
of argument we will take to be negative, distributed evenly throughout its 
volume. The body is considered to be porous, so that both multivalent and 
monovalent ions can reach and pass through its bulk. Finally, we assume that 
the solution occupies a volume much larger than that of the sphere, so that 
any perturbations in the electric potential or ion concentration caused by the 
sphere go to zero at large distances. A schematic of the model is given in 
Fig. [TJ The choice of a spherical geometry is arbitrary, but we note that for a 
cylindrical charged body the form of the solution is the same. The cylindrical 
case will be discussed more completely at the end of the following subsection. 
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It has been shown through numerical experiments [10] and physical arguments 
[TT] that monovalent ions respond linearly to changes in electric potential and 
therefore do not take part in the formation of strongly-correlated liquids that 
produces charge inversion. On the other hand, the multivalent counterions 
(of valency +Z, hereafter "Z-ions") may become strongly correlated to each 
other as they adsorb into the volume of the sphere. We assume that for small 
deviations from the neutral point, adsorbed Z-ions are bound to the sphere by 
the correlation energy /x, which we take to be constant for a given bare charge 
density and ion valency Z . The limitations of this assumption are examined in 
the Discussion section, but for small deviations of the counterion concentration 
from the neutral point our assumption is justified, as the properties of the 
strongly correlated liquid do not change significantly. 

Finally, we assume that the concentration of monovalents and multivalent 
salts are fixed and constant far from the sphere. Motivated by biological sit- 
uations, we will consider the case where the monovalent salt is much more 
prevalent than multivalent salt, i.e. the monovalent salt concentration iVi )00 
is much greater than the Z-ion concentration Nz )0 o- This assumption leads to 
the simplification that the Debye-Hiickel screening radius 



is determined solely by the concentration of monovalent salt, N = Ni j0a . Here, 
e ~ 80e is the permittivity of water, ksF is the thermal energy, and e is the 
proton charge. In biological situations, where salt concentration is typically on 
the order of 100 mM, the screening radius r s is less than 1 nm. So our analysis 
assumes that r s -C R. 



2.2 Free energy description 

When the sphere has a net charge, there is an expansive pressure acting on 
its surface, and this pressure can be computed by first putting together a full 
description of the Helmholtz free energy F of the system. Specifically, the 
pressure P at a given Z-ion concentration is 



where V oc R 3 is the volume of the sphere. We are interested in deviations in 
Z-ion concentration away from the neutral point, defined as the concentration 
at which the net charge on the sphere is zero. It is therefore natural to 
define this neutral point as our reference for free energy. 






(2) 
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In order to determine the free energy, it is necessary to solve for the electric 
potential in the vicinity of the sphere. We begin by assuming that a certain 
net charge Q = Q + MZe exists within the sphere's volume, where Q is the 
(negative) bare charge and M is the number of adsorbed Z-ions. In actuality, 
the existence of a nonzero net charge represents a deviation in the number of 
adsorbed Z-ions from the neutral point Mq = \Qo\/Ze, and the free energy 
associated with this deviation will be addressed momentarily. For now, how- 
ever, we assume only that the sphere has some net charge Q, and we use it to 
derive the electric potential. 

In the limit of large screening, r s <C R, the potential <fi relative to infinity is 
nearly constant within the volume of the sphere and can be approximated as 



where p is the volume density of the net charge, p = Q/^nR 3 . This result 
allows us to define a an effective capacitance of the sphere 

_ AireR 3 

If e(j)/kBT < 1 so that the linearized Poisson-Boltzmann equation may be 
used, the free energy F so i of the ionic solution is well-known [T2lfT3] : 

Fsoi = l -Q(t>- (5) 



Until now we have ignored the free energy associated with the Z-ions adsorbing 
to the sphere, taking for granted some net fixed charge Q. We now consider 
the process of adsorption/desorption that leads to the net charge. There is 
a change in both energy and entropy corresponding to the movement of Z- 
ions from solution to their adsorbed state within the sphere's volume. Let us 
therefore ascribe some free energy F a d s to this process so that the total free 
energy 

F = F sol + F ads . (6) 

A charge Q indicates that a total number M — M = Q/Ze Z-ions have 
adsorbed (or desorbed) to the sphere relative to the neutral point. If this 
number is small compared to the total number adsorbed onto the sphere, then 
we can consider the correlation energy p and bulk counterion density Nz,buik 
to be constant. The total free energy of adsorption is then 

Fads = -tt-H + -^-k B T\n (N z ,buik/N z ,oo) , (7) 
Ze Ze 

where the first term on the right-hand side represents the energy gain of the 
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Z-ions and the second term indicates the entropic contribution resulting from 
Z-ions being confined within the sphere at a density N Z) buik- Eq. ((Tj) can be 
rearranged to read 

F ^ = .Q k ^J N ^v(-Mlk B T) y 

Ze \ N z ,oo J 

Previous authors [7] have noted that the quantity Nz,buik exp(— \h\/}zbT) = 
Nz,o constitutes a new boundary condition for the Poisson-Boltzmann equa- 
tion in problems where multivalent ions condense onto a charged surface. Nz,o 
is a sort of "vapor concentration", giving the value of counterion concentra- 
tion just above a charged surface with correlation energy \i. In our case, this 
value can be used to define a balance between N ZyOQ and the concentration of 
adsorbed Z-ions using the Boltzmann distribution: 

N Zfi = N Zt00 exp {-Ze(j>/k B T) (9) 

or 



Eq. ( TTUj) suggests that the electric potential inside the sphere is a direct result 
of the concentration gradient between the Z-ions in solution N Z)OC and the 
"vapor concentration" N z<0 . The neutral point must therefore occur when 
N z ,oo = N z ,o- This result allows us to conclude that the electric potential and 
the amount of overcharge (by Eq. (JBJ)) are driven by the ratio of the Z-ion 
concentration to the neutral point concentration. 



Substituting Eq. and (1T0]) into Eq. flS]), we obtain the total free energy 
from adsorption as simply 

F ads = ~Q<P, (11) 

so that the free energy of the entire system is just 

F = F sol + F ads = - l -Q(f ) (12) 

or 

F=-- ^. (13) 

It is worth noting that Eq. (j!2p is identical in form to the free energy of 
a capacitor held at constant potential by a voltage source. In this way our 
sphere is like a capacitor (as in Eq. (TjJ ) held at constant potential by a 
concentration ratio Nz,ao/N z of Z-ions. 
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Using Eq. (j2J), we can solve for the pressure and conclude that 



P 



( 



dF 
dV 



) 



N z , 



dF 
8V 




(14) 



In other words, the electrostatic potential of the sphere remains constant as 
it expands, leading to a positive (expansive) pressure that is proportional to 
cf) 2 . In terms of counterion concentration, 



so the sphere experiences a parabolically-increasing expansive pressure when- 
ever the Z-ion concentration deviates from the neutral point. We note that the 
pressure in Eq. ( Fl4|) can be interpreted as the energy density of the electric 
field at the surface of the sphere. 

We now briefly consider the case where, rather than assuming a spherical 
geometry, we take the charged body to be a cylinder of length L whose radius 
a <C L is free to expand as multivalent counterions bind to its volume. We can 
solve for the potential using the method of Eq. (J3J) to get 4>cyi — Q/enL-r 2 /a 2 . 
The pressure is therefore 



That is, the outward pressure for a cylinder has the same form as for a sphere. 
This result is not particularly surprising, since for small r s the potential is con- 
stant throughout the body and therefore the net charge Q and the monovalent 
salt are distributed isotropically throughout its interior. In this sense the so- 
lution of charges within the body behaves like a Pascal fluid, and the outward 
pressure it creates (Eqs. (fill) and (fIB]) ) should be uniform and independent of 
the shape of the object containing it. 

2.3 Elastic expansion 

If some degree of elasticity is ascribed to the sphere (or cylinder), then the 
outward pressure given by Eq. (Tl5]) causes it to swell. For the DNA bundles 
described in the Introduction, we may expect elastic behavior as the result of 
a balance between short-range forces within the bundle: correlation-induced 
attraction and steric repulsion. When the bundle is completely neutral, its 
equilibrium interhelical spacing corresponds to a minimum in the total free 
energy. We may expand the free energy by Taylor series to second order in the 
interhelical spacing and arrive at a potential that grows parabolically as the 




(15) 




(16) 
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spacing deviates from this minimum (its value at the neutral point). When 
the bundle becomes under- or over- charged by adsorption/desorption of coun- 
terions, the existence of a net charge creates an electrostatic pressure which 
produces radial swelling in this parabolic elastic potential. For polyelectrolyte 
brushes, it is the brush height that swells in an elastic potential. When the 
brush is neutralized, it achieves a minimum height that corresponds to a min- 
imum of the total free energy (which includes electrostatic energy as well as 
conformational entropy). The free energy near this minimum may also be ex- 
panded by Taylor series to second order in the brush height, resulting in a 
parabolic potential. 



Our assumption of a parabolic potential for small deviations from the neu- 
tral point allows us to write down the free energy Fb associated with the 
compressibility of the body as 



B 



(17) 



where B is the bulk modulus and x is the fractional change in volume, x = 
AV/Vo- If we add this term to the free energy for our sphere we can calculate 
its swollen size using the equilibrium condition 



0. 



Let the swollen radius R be written as R = Rq + AR, where Rq is the radius 
at the neutral point. If AR is small compared to R, so that the change in 
volume AV oc AR, Eq. ( TP81 can be evaluated as 




AR 

Rn 



Vq e ,2 / 2 

6 /rz. 



(19) 



In other words, the swollen radius is directly proportional to the outward 
pressure (Eqs. (TH1) and ffT5]) ): 



AR oc 



In 



' N 



1 2 



Z,oo 



N. 



(0) 



(20) 
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3 Discussion 



3.1 Comparison with existing data 

Fig. [2] shows a typical plot of swollen radius versus counterion concentration. 
Overlaid on the plot are data from an experiment in Ref. [5J, where the in- 
terhelical spacing of closely-packed DNA condensates was measured via x-ray 
diffraction for various concentrations of spermine salt (valency 4+). In this 
experiment, the addition of spermine to a solution of short DNA fragments 
was seen to induce the formation of dense, ordered DNA condensates, as de- 
picted in the inset of Fig. [2J The interaxial spacing a between parallel DNA 
"rods", which in this experiment constitutes the swollen dimension, was ex- 
tracted via x-ray diffraction over a range of spermine salt concentration. The 
data is plotted in Fig. [2] as the ratio a/a^°\ where is the observed minimum 
spacing. 




Fig. 2. Swollen size of a charged body in a solution of multivalent counterions. 
Experimental data, shown as solid dots, is taken from Ref. [6], where adsorption 
of spermine salt (4+) produces a non-monotonic change in the interhelical spacing 
of condensed DNA aggregates (shown schematically in the inset, reprinted with 
permission from [6] ) . The value of is taken from Ref. [To] . The elastic constant 
as defined in Eq. (|17p is chosen to match the experimental data, and it is similar to 
the elastic constant predicted by Ref. [T7] . 

The value of the neutralizing concentration reported in this study (~ 3 
niM) is in strikingly good agreement with the predicted value for DNA- 
spermine condensates from Ref. [15] (3.2 mM), where reentrant condensation 
of DNA was described in the context of correlations between adsorbed multi- 
valent cations. The "elastic constant" of the parabola in Fig.[2]was determined 
by fitting the experimental data to the form of Eq. (1191) . The value of this con- 
stant compares favorably with other theories of correlation-induced attraction, 



10 



as described in the following section. 

Data for swelling of polyelectrolyte brushes, taken from Refs. [3j and [5] also 
shows non-monotonic dependence on salt conentration. Both studies reported 
measurements of the height of long polyacid brushes as divalent salts like 
Ca 2+ and Mg 2+ were added. The brushes were seen to collapse abruptly when 
low levels of the salts were added, and then to swell noticeably over a certain 
range of higher salt concentration. There is some disagreement between the two 
experiments on the degree of swelling in the brushes, and the ability of divalent 
ions to form a strongly-correlated liquid is marginal, but the convincingly 
non-monotonic dependence on ion concentration is encouraging. The authors 
of Ref. [3J also report measurements of the brush height as Al 3+ ions are 
added to solution. They find that, over the studied range of aluminum salt 
concentration, the brush is in a collapsed state but swells noticeably with 
increased amounts of Al 3+ . Because of the exponential dependence of on 
Z, we suspect that this range (0.1 to 100 mM of AICI3) falls above the neutral 
point, so that the net charge of the brush is already inverted and further 
addition of trivalent salt produces outward swelling. 



3.2 Limitations of elastic expansion 

In deriving the parabolic form of Eq. (|20|) . we assumed that the properties 
of the adsorbed strongly-correlated liquid were unchanged as the sphere ex- 
panded. The correlation energy /x was held constant, and the sphere was as- 
sumed to respond elastically as the radius increased. In reality, increasing the 
size of the sphere causes the adsorbed ions to become more distant, lowering 
the correlation energy between them. This effect leads to significant deviations 
from parabolic behavior when the Z-ion concentration is far from the neutral 
point. 

In the case of coiled DNA or other polyelectrolytes, the correlation energy is 
in fact what holds the polymer in its condensed phase. That is, our assumed 
"elastic force" in DNA is the result of a balance between steric repulsion and 
correlation- induced attraction among interior charges, and therefore depends 
directly on /j. Thus as the polymer swells and the correlation energy between 
adsorbed ions decreases, the contractile "bulk modulus" of the condensate 
must also decrease. In this case the polymer will swell faster than predicted 
by the elastic model. For many polyelectrolyte condensates there will be some 
degree of swelling at which the repulsive electrostatic energy becomes greater 
than the maximum attractive correlation energy, causing the polyelectrolyte to 
"blow up" . This phenomenon results in sharp transitions between condensed 
and dissolved phases for DNA in solution, as has been described theoretically 

USES]. 
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The way in which the "elasticity" of a charged body changes with its size 
may be complicated, and in general it depends on the physical structure of 
the body. However, we may gain insight into changes in the elastic energy Fb, 
and therefore of the limitations of parabolic swelling behavior, by considering a 
basic example of correlation-induced attraction. The authors of Ref. [17] have 
assembled a theoretical description of attraction between like-charged plates 
in a solution of counterions that may be relevant to our situation. Using Monte 
Carlo simulations, they find scaled numerical values for the attractive pressure 
between parallel plates as a function of their separation. They demonstrate 
that the pressure falls off quickly when the separation between plates grows to 
the order of the average distance between adsorbed ions. This effect is similar 
to what we might expect for a condensed polyelectrolyte like DNA: when the 
condensate swells to a certain characteristic size the correlation energy will 
drop off quickly, leading to more rapid swelling of the polyelectrolyte and 
eventual resolubilization. 



We can fit the numerical solution for attractive pressure between plates [UJ to 
the situation of a crystalline DNA condensate (as in Ref.JB]), where we expect 
the attractive force between parallel DNA helices to decrease as they are pulled 
apart by swelling. Using the data shown in Fig. [2j we can evaluate the swollen 
size of a DNA condensate whose elasticity follows this "parallel plate" model 
and predict when it might "blow up" due to swelling. We scale the results for 
pressure between plates in Ref. [T7j (their Fig. 3) so that the bulk modulus 
B = —VdP/dV at equilibrium matches the observed bulk modulus for DNA 
condensates near the neutral point. In this way we determine an equation for 
the bulk modulus of a DNA condensate bundle as a function of the radius, 
B = B(AR/ R ). A numerical iteration of Eqs. ([3]) and (fT9j) can then be 
performed at a given Z-ion concentration in order to find the equilibrium size. 



The result of this procedure is shown in Fig. [3j and highlights the sharp 
deviation from parabolic behavior when the Z-ion concentration is far from 
the neutral point. It is worth noting that the bulk modulus resulting from 
fitting experimental data (6.5 MPa) is similar to what one might estimate 
[T8] based on the parallel-plate model of Ref.jTT] for the equilibrium point (13 
MPa). The experimental data presented in Fig. [21 however, does not deviate 
from the elastic model as strongly as the prediction of Fig. [3], and this difference 
may be attributed to the increased flexibility and conformational degrees of 
freedom of the DNA segments as compared to two parallel plates. In fact, the 
dotted line in Figj3] represents a sort of "worst-case scenario" for deviation 
from the simple parabolic description of Eq. (|20|) . where the system has no 
freedom to change conformation and the correlation energy must decrease 
abruptly with swollen size. 
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Fig. 3. Swelling as a function of Z-ion concentration, using an elastic model (solid 
line) and a "similarly charged plate" model (dashed line) where elastic force de- 
creases with swollen size. The bulk modulus is taken from the data of Ref. [6] as 
in Fig. and the similarly charged plate model is adapted and scaled from Ref. 
[17j . The charged plate model deviates significantly from the elastic model when 
the Z-ion concentration is far from the neutral point, driven by a decrease in the 
correlation energy between adsorbed ions. 

4 Concluding remarks 

We have shown in this paper how a porous macroion may exhibit non-monotonic 
swelling behavior as multivalent ions are added to the surrounding solution. 
We have presented a simple model of swelling near the neutral point based on 
the concept of correlation-induced adsorption of ions, and concluded that the 
swollen size of the macroion grows as the square of the logarithm of the ion 
concentration. 

As a final comment, we note that the applicability of our conclusion depends 
strongly on the ability of adsorbed multivalent counterions to form a strongly- 
correlated liquid. The degree to which this is possible depends in general on 
the background to which the multivalent counterions can adsorb. If the "bare 
charges" distributed throughout the macromolecule are too distant, then such 
a strongly-correlated liquid will not be possible; the distance between adsorbed 
Z-ions must remain smaller than the screening radius r s or there will be no 
significant correlation energy. The valency of the counterions also plays an im- 
portant role: the correlation energy /x of condensed counterions grows directly 
[7j with Z, and the concentration of multivalent ions required to produce over- 
charging decreases exponentially with /j. Theoretical estimates [7f20] indicate 
that our approach is valid for Z > 3, while the case of Z = 2 is marginal. A 
recent paper [19] has specifically considered the case of divalent ions adsorb- 
ing to a flexible polyelectrolyte, including the effects of correlations between 
adsorbed ions along the polyelectrolyte backbone. The authors conclude that 
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even a modest concentration of divalent salt can produce overcharging and 
non-monotonic swelling at room temperature. 
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List of Figures 



1 Schematic of the system to be studied. A porous, spherical, 
negatively charged body sits in a solution of positive 
multivalent ions (denoted by +Z's) and dissolved monovalent 
salt (+'s and -'s). Monovalent salt passes freely through the 
sphere while the Z-ions may adsorb to its volume, forming a 
strongly-correlated liquid. 

2 Swollen size of a charged body in a solution of multivalent 
counterions. Experimental data, shown as solid dots, is 
taken from Ref. [6], where adsorption of spermine salt (4+) 
produces a non-monotonic change in the interhelical spacing of 
condensed DNA aggregates (shown schematically in the inset, 
reprinted with permission from [6]). The value of is taken 
from Ref. [15]. The elastic constant as defined in Eq. (ITT)) is 
chosen to match the experimental data, and it is similar to the 
elastic constant predicted by Ref. [TT| . 

3 Swelling as a function of Z-ion concentration, using an elastic 
model (solid line) and a "similarly charged plate" model 
(dashed line) where elastic force decreases with swollen size. 
The bulk modulus is taken from the data of Ref. [6] as in Fig. 
El and the similarly charged plate model is adapted and scaled 
from Ref. [TTJ. The charged plate model deviates significantly 
from the elastic model when the Z-ion concentration is far 
from the neutral point, driven by a decrease in the correlation 
energy between adsorbed ions. 
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